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ABSTRACT 

The fraction of ionizing photons that escape from young star-forming galaxies is one of the largest 
uncertainties in determining the role of galaxies in cosmic reionization. Yet traditional techniques for 
measuring this fraction are inapplicable at the redshifts of interest due to foreground screening by the 
Lyman a forest. In an earlier study, we demonstrated a reduction in the equivalent width of low-ionization 
absorption lines in composite spectra of Lyman break galaxies at z ~ 4 compared to similar measures at 
z ~ 3. This might imply a lower covering fraction of neutral gas and hence an increase with redshift in 
the escape fraction of ionizing photons. However, our spectral resolution was inadequate to differentiate 
between several alternative explanations, including changes with redshift in the outflow kinematics. Here 
we present higher quality spectra of 3 gravitationally lensed Lyman break galaxies at z ~ 4 with a spectral 
resolution sufficient to break this degeneracy of interpretation. We present a method for deriving the 
covering fraction of low-ionization gas as a function of outflow velocity and compare the results with similar 
quality data taken for galaxies at lower redshift. We find a significant trend of lower covering fractions 
of low-ionization gas for galaxies with strong Lya emission. In combination with the demographic trends 
of Lya emission with redshift from our earlier work, our results provide new evidence for a reduction in 
the average H I covering fraction, and hence an increase in the escape fraction of ionizing radiation from 
Lyman break galaxies, with redshift. 

Subject headings: cosmology: reionization — galaxies: evolution — galaxies: formation — galaxies: ISM 



1. INTRODUCTION 

Star forming galaxies are the leading candidate for the 
source of ultraviolet photons required to reionize the uni- 
verse. Several lines of evidence indicate that reionization 
was underway b y z = 1 1 and end ed a few hundred My r later 
at z ~ 6-7 fe.g-ISchenker et al.ll2012tlMortlock et al.|[2011t 
iHinshaw et al.ll2012l ). Deep near-IR imaging with the Hub- 
ble Space Telescope has provided good constraints on the 
UV luminosity density of star forming galaxies during the 
reionization epoch (e.g.. lEllis et al.ll2~013HOesch et al.ll2013[ ) 
which indicate that the ongoing star formation is likely ca- 
pable of producing the required ionizing flux. However, it is 
unclear whether this radiation is able to escape from galax- 
ies and actually ionize the intergalactic medium (IGM). 
Based on estimates of the total UV luminosity density and 
IGM clumping factor, the required escap e fraction f PRn of 
hydro gen- ionizing photons is > 0.2 (e.g., iRobertson et al.1 
2013). The precise value of f esc is a key uncertainty in 
determining the role of galaxies in reionization. 

Direct measurements of f esc during the reionization epoch 
are essentially impossible, not only because of the faint ap- 
parent luminosity, but also because foreground IGM at- 
tenuates the ionizing flux to undetectable levels even at 
z > 4. Direct imaging and composite spectra of galaxies at 
lower redshift has established a modest aver age f P . sr ~ 0.05 
for L yman break galaxies (LBGs) at z = 3 fBog osavljevid 
2010). If star forming galaxies dominated reionization, f esc 
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must have been higher at earlier times. This paper is con- 
cerned with improving the constraints on f esc at higher 
redshifts for which direct measurements are not practical. 
Our methodology is as follows: f esc is set by the areal cov- 
ering fraction of hot stars by H I such that f esc = 1 — f c , 
and f c can be inferred from intermediate dispersion spec- 
troscopy of interstellar UV absorption lines. The difficulty, 
of course, is that LBGs at redshift z > 3 are very faint, 
so securing suitably high quality absorption line spectra of 
individual examp les is a very chall enging proposition. In 
an earlier paper (| Jones et al.l 120121) we therefore analyzed 
the average properties of LBGs at z ~ 4-5 derived from 
composit e spectra in a manne r similar to that pioneered at 
2-3 bv lShaplev etUI (pOOl . 

A particular motivation for the present study was the dis- 
covery of a marked reduction with increasing redshift in the 
equivalent width of low-ionization absorption lines at fixed 
UV luminosity, suggestive of changes in either t he kinematic 
profil e or the covering fraction of neutral gas (jJones et al.1 
l2012f l. Additionally, spectra of individual galaxies show 
a trend of increasing Lya emission equivalent width with 
redshift which we arg ued could reflect evolution in the 
H I covering fraction (Stark et al.ll2010| ), especially given 
the strong correlation of Lya and low-ionizatio n absorp- 
tion lines (jShaplev et al.l 120031: [Tones et al.ll2012D . A lower 
covering fraction for the neutral gas within typical LBGs 
would be particularly important as it could imply a higher 
escape fracti on of ionizing phot ons. The composite spectra 
discussed bv lJones et al.l ()2012f ) did not have adequate reso- 
lution to distinguish the effects of reduced covering fraction 
and kinematics, and so the present paper takes this investi- 
gation one step further by attempting to resolve this impor- 
tant ambiguity. Here we present higher resolution spectra 
of 3 gravitationally lensed z ~ 4 LBGs. Although their un- 



lensed luminos ities are typ i cal of the constituent galaxies 
comprising the lJones et al.l (|2012f) composite, their individ- 
ual lensed magnitudes are much brighter enabling compa- 
rable signal to noise to the stack of LBGs discussed in our 
earlier paper. 

Throughout the paper we adopt a flat ACDM cosmology 
with fi A = 0.7, VL M = 0.3, and H = 70 kms^ Mpc' 1 . All 
magnitudes are in the AB system (|Okelll974l ). 

2. GRAVITATIONALLY-LENSED Z ~ 4 GALAXIES 

The spectr a com prising the composite published by 
Uones et al.1 (J2012I ) were taken with the 600 line mm -1 
DEIMOS grating with a resolution of ~3.5 A, although un- 
certain systemic redshifts led to reduced resolution of the 
composite spectrum (corresponding to a velocity resolution 
of ~ 450 kms -1 FWHM). The composite comprised galax- 
ies with z' AB = 24-26 with 90% completeness to z' AB = 25 
at z ~ 4. The stacked spectrum reveals multiple low ioniza- 
tion lines such as Si II A1260, O I A1302 + Si II A1304, and 
C II A1334 in t he region where th e signal/noise is optimal. 
As discussed bv lJones et al.1 <|2012f ). these lines are normally 
saturated so the line depth at a given velocity provides a 
measure of the areal covering fraction f c of O and B stars 
by neutral H I gas along the line of sight. 

Typical LBGs are too faint for detailed line profile stud- 
ies even at z — 2, but strong gravitational lensing can 
boost the brightness of representative examples making 
such studies of individual sources a practical proposition. 
Studies of several lensed z ~ 2-3 LBGs have found that ab- 
sorption velocities of low-ionization metal transitions range 
from ~ —1000 to +500 kms" 1 with typical line centroids 
y ~ -200 kms" 1 (IPettini et al|[200l IQuider etliHl2T)09l 
120101 iDessauges-Zavadsky et al.H2010D. The mean lo w ion- 
ization absorption velocity in the I Jones et al.l ()2012f ) com- 
posite is similar, vlis = —190 kms -1 . 

In a similar fashion, for the present analysis we have 
located 3 gravitationally-lensed LBGs at z ~ 4. Two 
are independent sources lensed by the well-studied clus- 
ter Abell 2390 and the third was located in the cluster 
J1621+0607 in the Sloan Digital Sky Survey. A2390_H3 
and H5 represent two distinct highly-elongated pairs of 
lensed images that were spectroscopic ally confirmed to be 
at different r edshifts by , resp ectively, iFrve fc Broadhurstl 
(J1998D and iPello et~aTT (|1999T) . The tangential arc sys- 
tem in J162 1 +060 7 was spectroscopically confirmed by 
IBavliss et al.l (|2011f) . T he gravitational m agnification fac- 
tor is ~10 in each case (jPello et al.lll999|) . We summarize 
the key properties in Table Q] The redshifts and absolute 
UV lumino s ities a re representative of sources studied by 
Uones et~a l. (2012). In addition, for comparison purposes, 
we include an analysis of high quality spectra for 3 further 
lensed z = 2-3 sources (courtesy of M. Pettini) in Table CD 
These include t he 'Horseshoe ' (z=2.38, IQuider et al] l2010). 
cB58 (z =2.73. IPettini et all 120021 ) and the 'Cosmic Eye' 
(2=3.07, IQuider et al.H2009D . 

Spectra were taken with the 1200 line mm -1 DEIMOS grat- 
ing during two runs in October 2011 and June 2012. This 
provides a resolution of ~1.7 A corresponding to a veloc- 
ity resolution of ~70 km s -1 , considerably be tter than for 
the composite discussed by Uones et al.l (|2012l ). Spectra of 
each galaxy covered wavelengths corresponding to at least 
1175—1675 A in the rest frame. The lensed sources in Abell 



2390 were observed simultaneously with a multi-slit mask 
that sampled two images of each source. Seeing varied be- 
tween Of! 4 and l'/4 FWHM during the observations, and 
the bulk of the data used has seeing in the range 0'.'7-0'.'9. 
Some exposures (~ 10%) were affected by cirrus and are 
not included in the final addition. Total observing times 
for the final spectra are given in Table [TJ 
The DEIMOS spectra were reduced and calibrated using 
the Spe c2D pipeline follow ing the techniques discussed in 
detail bv lStark etaP (j2010f l. In the case of A2390JI3, care 
was taken to ensure that the extracted spectrum was not 
contaminated by light from a nearby cluster member. Data 
from the October 2011 and June 2012 runs were reduced 
separately and the resulting one-dimensional spectra were 
combined with an inverse-variance weighted mean. Spec- 
tra of J1621 are affected by poor sky subtraction residuals, 
while the Abell 2390 arc spectra are of excellent quality. 
Spectra of different images of the Abell 2390 arcs were 
scaled to the same flux level before combining to a com- 
mon wavelength scale with 0.7 A pixels, roughly Nyquist 
sampled. The final spectra, shown in Figure [TJ reach an 
average continuum S/N per 70 kms -1 resolution element 
of 5 for J1621, 9 for A2390.H3, and 10 for A2390_H5 over 
the rest-frame wavelength range 1250-1650 A. This is com- 
parab le to that in the composite spectrum in Uones et al.l 
(2012) which has S/N equivalent to ~ 10 at the improved 
resolution of 70 kms^of our new data. 

3. ANALYSIS 

3.1. Systemic Redshift 

Accurate systemic redshifts are required in order to ex- 
amine the kinematics of gas seen in absorption and the 
techniques fo r esti mating these are discussed in detail in 
Uones et al.l (|2012t ). This is straightforward when neb- 
ular emission lines are visible, such as is the case in 
both J1621+0607 (O in] AA1661,6) and A2390JH5 (O in] 
AA1661,6, He n A1640, C iv AA1548,51). The strong emis- 
sion from highly ionized species such as He II and C iv 
seen in A2390_H5 is uncommon but has been observed 
in some high redshift starburst galaxies and signifies an 
extre mely young, metal-poor, and hot st ellar population 
(e.g. iFosburvet all 120031: lErb et al.ll2010D . Alternatively 
they may signify the presence of an active galactic nucleus, 
but the narrow line widths (50 — 185 kms -1 FWHM, cor- 
rected for instrumental resolution) suggest an origin in star- 
forming Hn regions. 

No appropriate features are detected in the spectrum 
of A2390_H3 and so we estimate the systemic redshift 
from low-iqnizatio n absorption lines using the method of 
Uones et alJlpOll ). This gives z = 4.043±0.002 = z /s + 190 
kms -1 , with uncertainty dominated by an rms difference 
~ 125 kms -1 between redshifts derived fro m absorption 
lines and that obtained from nebular emission (jSteidel et al.l 
2010). The adopted systemic arc redshifts are listed in Ta- 
ble [Q 

3.2. Low- Ionization Covering Fraction 

Ideally we would measure the covering fraction of neu- 
tral hydrogen directly from spatially resolved H I absorp- 
tion. However, the only available transition (Lya) is domi- 
nated by strong emission with net equivalent width W\, ya = 
20 — 100 A, and the observed line profile is complicated by 
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resona nt scattering in the extended circumgalactic medium 
(CGM; ISteidel et al.ll201ll ). These effects are apparent from 
the Lya line profiles which show redshifted emission as well 
as strong absorption arising from both the CGM and Lya 
forest (Figure [T]) . We therefore estimate the covering frac- 
tion of neutral hydrogen from absorption lines of heavier 
low-ionization species which arise in H I gas, i.e., those with 
ionization potentials less than 1 Rydberg. 

The covering fraction of any ion is related to its absorption 
line optical depth r and residual intensity I via 



T = 1 " /c(l - O 



(1) 



where Iq is the continuum level. Optical depth is in turn 
related to column density as 



m e c 



/A 
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3.768 x 10 14 



(2) 



where / is the ion oscillator strength, A is the transition 
wavelength expressed in A, and N is the ion column density 
in cm~ 2 (kms -1 ) -1 . Combining equations [T] and [2] yields 
an expression for f c as a function of / and N. In cases 
where two or more transitions are measured for the same 
ion, from the same ground state, with different values of 
/A, it is possible to solve these equations for N and f c . In 
the following analysis we will treat all variables as functions 
of velocity, i.e., f c (v). 

For the low ionization species of interest, our spectra cover 
three such transitions of Si n at 1260, 1304, and 1526 A 
which we use to measure the covering fraction as a func- 
tion of velocity f c (v), for each galaxy. Si n A1304 is only 
used in the velocity range v > —200 kms -1 where it is not 
contaminated by O I A 1302. We resample the spectrum 
of each transition to a common velocity scale, and find the 
values of N and f c which minimize the least-square residual 
X 2 = 52(Iobs - In J J 2 /crl bs in each velocity bin. We addi- 
tionally find the range of N and f c for which x 2 is within 1 
of the minimum value, and adopt this as the la uncertainty. 

The best-fit f c and uncertainty calculated for each arc are 
shown as a function of velocity in Figure [2] Since Si n 
is the dominant ion of silicon in H I gas, this is approxi- 
mately equal to the covering fraction of H I (provided that 
it is enriched with Si) which impedes the escape of ionizing 
radiation. 

3.3. Average Low- Ionization Absorption Profile 

A simple and complementary alternative to the method 
outlined in Section 13.21 is to estimate the covering fraction 
from saturated transitions. In cases where r» 1, Equa- 
tion [T] simplifies to 



u = i - i/h. 



(3) 



Several of the strongest absorption lines covered by our 
spectra are typically saturated, including: Si II A1260, O I 
A1302, Si ii A1304, C n A1334, and Si n A1526 which are 
all tracers of H I gas. In order to minimize the statistical 
uncertainty we calculate the average intensity of these tran- 
sitions as a function of velocity using an inverse-variance 
weighted mean, taking care not to use the wavelength re- 
gion where O I A1302 and Si n A1304 are blended (roughly 
—300 < v < —200 kms -1 depending on the kinematics 



of each source). These profiles are shown in Figure [5] to- 
gether with the covering fraction measured from Si II. We 
note that the covering fraction derived from Equation [3] is 
a strict lower limit. 

4. RESULTS 

This work was motivated in large part by the need to disen- 
tangle kinematics and covering fractions of absorbing gas. 
In particular, we seek to explain the extent to which de- 
creased absorption line equivalent widt hs measured from 
composite spectra in our earlier work ([Jones et al.M2012T) 
result from changes in gas kinematics compared to cover- 
ing fractions, and the implications of this result for the es- 
cape of ionizing radiation. We are limited in examining the 
redshift evolution of these properties by the small number 
of sources with suitable spectra, and we caution that this 
sample is not necessarily representative of the LBG popula- 
tion at these redshifts. Nonetheless we can examine general 
trends within the existing d ata from this work and others in 
the literature dPettini et al. 2002; Quider et al.ll2009t 120101: 
Dcssaugc s-Zavadskv et al.ll2010( ) . Following the methods in 
Section [3] we show the average absorption profiles and Si II 
covering fractions (derived from unblended transitions at 
1260, 1526, and 1808 A) of well-studied z = 2-3 galax- 
ies for comparison in Figure [2j To quantify trends in the 
absorption line profiles, the velocity extent and maximum 
absorption depth for each galaxy are shown as a function 
of redshift in Figure [U 

In the case of Abell 2390, we take the most conservative ap- 
proach noting that the interpretation of the absorbing gas 
is complicated by the physical proximity and similar red- 
shift of the t wo ar cs. Their projected separation is ~ 70 kpc 
(jPello et al.lll999t ) and it is uncle ar which sour c e lies in the 
foreground. At lower redshifts, iSteidel et al.l (j2010f ) have 
shown that low-ionization absorption seen in a background 
source at b — 70 kpc has a detectable average equivalent 
width ~ 0.4 A for the transitions of interest. Since we lack 
specific information about the 3-D geometry, the following 
analysis does not include any contribution from this effect. 
If anything, our results will overestimate the true covering 
fraction and therefore yield a more conservative constraint 
on the escape fraction. 

4.1. Kinematics 

The kinematics of foreground low-ionization gas are re- 
vealed in the absorption line profiles shown in Figure [5J 
In all cases we see significant blueshifted absorption indi- 
cating outflows, a s expected give n the high star formation 
surface densities (jHeckma nl I2002D . The maximum outflow 
velocity at which absorption is detected is —700 kms -1 in 
A2390_H3, with an uncertainty of ~ 125 kms -1 since we 
do not directly measure the systemic redshift. A2390_H5 
reveals weak absorption extending to —600 kms , seen 
also in higher-ionization Si iv and C iv lines, although it is 
only marginally detected at < —300 kms" 1 . The maximum 
outflow velocity in J1621 is —300 kms -1 . 

The outflowing low-ionization gas attains a somewhat lower 
(~ 30% on average) maximum velocity at higher redshift. 
This trend is not due to lower quality data as it remains ev- 
ident in Figure[5]if we consider alternative measures such as 
the FWHM or an absorption threshold at 25% of the con- 
tinuum flux. However all galaxies except J1621 have similar 



maximum velocities ranging from 600 — 800 kms -1 . Like- 
wise the extent of redshifted absorption is approximately 
+200 kins -1 for all sources, with the nota ble exception 
of the Cosmic Eye as discussed in detail by iQuider et al.1 
(|2010() : this indicates little difference in line broadening 
from rotation or other internal kinematic structure. There- 
fore, in this limited sample, the low-ionization gas kine- 
matics are similar with a somewhat lower average velocity 
extent at higher redshift. 

4.2. Covering Fraction 

The covering fraction of each galaxy as a function of gas 
velocity is estimated from the methods described in Sec- 
tion [3] and shown in Figure [2l Both methods are generally 
in good agreement indicating that Equation [3] is a valid 
approximation. Si II covering fractions derived for the Cos- 
mic Eye are systematically higher than indicated by the 
average absorption profile; this is largely an artifact caused 
by additional absorption at the w avelength of Si II A 1260 
from intervening gas at z = 2.66 (Qui der et al . 1 120101 ) . It 
is also apparent from Figure [2] that the Si II covering frac- 
tion is poorly constrained in regions of weak absorption 
due to the marginal significance of individual absorption 
lines. This is most problematic in the high-velocity wings. 
The strong anticorrelation between absorption line strength 
and covering fraction, a s well as results at lower redshift 
(jMartin fc Bouche1l2009f ). suggest that the most likely so- 
lution for such ambiguous cases is a low covering fraction 
of optically thick gas. We therefore opt to compare galax- 
ies on the basis of their average absorption line profile as 
this quantity is simpler to define and less susceptible to the 
uncertainties described above. Nonetheless the Si n results 
are an important verification that the average profile accu- 
rately traces f c . 

We can now compare the covering fractions measured at 
z = 4 with sources at lower redshift. Figure [3] shows the 
maximum absorption depth for each galaxy as a function 
of redshift. The z = 4 galaxies have maximum absorption 
depths corresponding to f c = 0.3-0.9, in each case occur- 
ring at v ~ —100 kms" 1 . There is a large scatter in Fig- 
ure |3] with o~(f Ct max) = 0.26 and no strong redshift depen- 
dence. Galaxies at z = 4 do, however, have covering frac- 
tions which are lower on average by 25% or Af Ctmax = 0.16 
compared to z = 2-3. 

4.3. Trends with Lya 

We now turn to trends with Lya equivalent width. Pre- 
vious sections focused on possible redshift evolution of 
low-ionization absorption lines in our quest to examine 
whether this may signify an increasing ionizing escape frac- 
tion. The connection with Lya is a natural one to ex- 
plore given there is a strong correlation betwee n its equiva- 
lent width WT.vry an d low-ionization absorption (jJones et al.1 
2012; iShapley et a l. 2003). Since our previous work has 
suggested that the distribution of W^v a for LBGs of a 
fixed luminosity increases w ith redshift (St ark et alj .2010. 
120111: iSchenker et al.l [20121. we can hope to derive infer- 
ences about the low ionization absorption in sources for 
which Lya measurements are now widely available. 
Each galaxy in Figure [3] is color-coded according to W^Lya- 
This value refers only to the equivalent width of Lya emis- 
sion, differing from the conventional net sum of emission 
and absorption. The maximum outflow velocity is lower 



on average in galaxies with stronger Lya emission, con- 
sistent with well-qua ntified results from composite spectra 
(|Shaplev et al J 120031 ) . More interestingly, Figure [3] reveals 
a trend of lower absorption depth (implying lower f c ) with 
stronger Lya emission at 3.5er significance. We show this 
relation in Figure [4] In contrast, the trend of lower av- 
erage f c at higher redshift has limited significance (l.ler) 
and is explained by Lya demographics within the sample. 
Although limited by the small sample size, this is an impor- 
tant first quantitative result at these redshifts. Since the 
frequency and equivalent width of Lya emission increases 
in LBGs at higher z = 3 — > 6, these results imply that 
the average covering fraction of low-ionization gas should 
decrease with redshift. 

4.4. Ionizing Escape Fraction 

Direct measurements of the ionizing flux are impractical at 
the redshifts of interest in this paper both because of the 
faint apparent magnitudes of LBGs and the high opacity 
of the Lya forest. Nonetheless we can provide important 
constraints on f esc using indirect tracers of H I. 
Before doing so, we consider the potential systematic un- 
certainties which may limit our ability to estimate the true 
value of f esc from metal absorption lines. Results at z ~ 3 
have shown that f esc is indeed dependent on low-ionization 
absorption strength (Bogosavlicvic 2010) although this re- 
lation is not one-to-one, likely due to the factors described 
below. While in general these preclude accurate estimates 
of the escape fraction, the maximum absorption depth (Fig- 
ures[2j[3]) is a valuable constraint on H I spatial homogeneity 
and sets a stringent upper limit on f esc . 

1. We measure covering fraction as a function of velocity, 
yet gas at different velocities may cover different spatial re- 
gions and we lack the spatial resolution needed to evaluate 
this effect. The maximum absorption depth at a given ve- 
locity is thus a lower limit on the total H I covering fraction 
and an upper limit on f esc . 

2. Metal- free H I is not detected. To date only two in- 
stances of metal-free ga s have been found at these redshifts 
(jFumagalli et al. 2011) and so this is likely insignificant. 
Again, this possibility implies that covering fractions mea- 
sured in Section [3] yield upper limits on f esc . 

3. Low column density gas will not be detected, although 
such gas will not affect f esc unless it has very low metallicity 
< 0.1 Zq. The optical depth of metal transitions used here 
compared to Lyman continuum at ~ 900 A is - J — = 1.7 for 

the weakest line (Si II A1304) and 10-20 for the strongest 
transitions (Si n A 1260, C II A1334, O I A1302) for solar 
abundance ratios ([Asplund et al.l [2009) . Galaxies in the 
z = 2-3 sample with measured interstellar abundance ratios 
have Z > 0.4 Z Q for the relevant elements, such that the 
attenuation of weaker Si lines is roughly equal to that of 
ionizing radiation. 

4. There may be narrow components with uniform cover- 
ing fraction which are spectroscopically unresolved. How- 
ever, such smooth absorption and column density profiles 
would require a remarkably regular spacing of discrete nar- 
row component s, and s o a partial c overing fraction app ears 
more likely (see lPettini et al.ll2002t IQuider et al.ll2009f ). 

5. Si II and C II are present in both H I and H n regions. 
We have therefore confirmed that O I, with essentially the 
same ionization potential as H I, gives consistent results in 
all cases. 
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6. We measure a covering fraction at 1260-1526 A, but 
the stars which are bright at these wavelengths do not 
necessarily emit at <912 A. Our measurements therefore 
correspond to a constraint on the relative escape fraction 
f esc = iLyc/^1500 as commonly used in the literature. 
However, the spatial distribution of £1500 is similar to ioniz- 
ing e mission as traced b y Balmer lines in z ~ 1 — 3 galaxies 
(e.g.. Uones et alJl2010t ). indicating that this effect may be 
minimal. 

With these caveats in mind, we now summarize what can 
be learned about the ionizing escape fraction at z = 4. The 
depth of low-ionization absorption lines gives upper limits 
f esc < 0.7 for A2390_H5, < 0.3 for A2390.H3, and < 0.1 for 
J1621. The true values are likely well below the upper limits 
and in thre e lower- redshi f t gala xies we can verify that this 
is the case. lOuider et al.l (|2010f ) show that the Cosmic Eye 
has an H 1 covering fraction of ~ 95% implying f esc < 0.05 
based on damped Lya absorption, more stringent than the 
f esc < 0.1 from our analysis. Secondly, no ionizing radi- 
ation is detected from the Horseshoe in deep UV imaging 
(B. Siana, private comm unication) nor from the spectrum of 
Q0000-D6 (f esc < 0.16; iGiallongo et al.l I2002T) despite low 
maximum covering fractions (Figures [2] [3]) . Non- uniform 
coverage of low-ionization metals is evidently a necessary, 
but not sufficient, condition for the escape of ionizing ra- 
diation. Results derived in this paper should therefore be 
strictly interpreted as upper limits. Nonetheless, our mea- 
surements at z — 4 are readily compatible with the value 
fear. ?", 0-2 required fo r galaxies to reionize the universe 
^Robertson et al.ll2013h . 



resolve this issue and provide better constraints on f c and 

J esc- 

We have successfully measured the kinematics and cover- 
ing fraction of low-ionization (H i) gas in z — 4 LBGs from 
high quality rest-UV spectra. Resulting upper limits on 
the ionizing escape fraction f esc < 1 — f c are readily consis- 
tent with that req uired for star-forming g alaxies to reionize 
the universe (e.g.. lRobertson et al.l l2013). Importantly, the 
new data enable these measurements at a time < 1 Gyr 
from the reionization epoch, much earlier than was previ- 
ously possible. We note that in order for galaxies to reion- 
ize the universe, the escape fraction must increase rapidly 
from {fesc) = 0.05 measured at z = 3 (iBogosavlievIcl feOlO) 
to a value f eS c ^ 0-2 at z = 7 (Roberts on et all 12013; 
iKuhlen fc Faucher-Giguerel 120121 ). While the trends with 
redshift are poorly constrained, the available data reveal 
a reduced covering fraction with increasing W^ yct indicat- 
ing that galaxies with moderate or strong Lya emission are 
likely to have larger fesc- This is supported by direct evi- 
dence at z = 3, where galaxies with detectable ionizing flux 
have stronger Lya emission and wea ker low-ionization ab - 
sorption than LBGs with lower f esc (jBogosavljevicI 12010) . 
Clearly the distribution of (metal-enriched) H I is patchier 
in galaxies with stronger Lya emission, and therefore more 
likely that ionizing radiation can escape. Since the fre- 
quency and strength of Lya emission in typi cal LBGs in- 
creases with redshift ([Stark et al.l l20ldL 1201 If ), our results 
provide new evidence that the covering fractions decrease 
(and therefore / esc increases) with redshift. 



5. DISCUSSION 

In our previous work ([Jones et al.ll2012T) we found a decrease 
in the low-ionization absorption equivalent width with red- 
shift in composite LBG spectra, but were unable to dis- 
tinguish whether gas kinematics and/or covering fraction 
were the cause. The new data presented in this paper now 
resolve this ambiguity. Although the sample of LBGs with 
high quality spectra at z > 2 is small and not necessar- 
ily representative, the present data show an approximately 
equal decrease of ~ 25% in the gas velocity and covering 
fraction with redshift. The much stronger dependence of 
low-ionization absorption with Lya emission appears to be 
predominantly due to the covering fraction of neutral hy- 
drogen. 

The main limitation of methods used in this work to con- 
strain f e sc is that we do not directly measure the fraction 
of ionizing emission covered by H I. While the spectrally- 
resolved covering fraction of heavy elements provides an 
important measure of the "patchiness" of H I, gas at dif- 
ferent velocities may cover different spatial regions and so 
we can derive only a lower (upper) limit on f c (/ esc ). Deep 
integral field spectroscopy with good spatial resolution may 
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b HST/WFPC2 F 


Sand et al.l (]2005l' 





c Derived from the spectrum at the mean wavelength of HST/WFPC2 F814W, A mcan = 8369.5 A. 

d For a magnification factor fi = 10. 

e This source is not gravitationally lensed. 



TABLE 2 
Spectroscopic Properties 
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a Emission component only. 

b Lower limit derived from absorption line depth. 
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On The Escape Fraction of Ionizing Photons in z ~4 Galaxies 
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Fig. 1. — Keck DEIMOS spectra of (top to bottom) A239CLH3, A 2390-H5, and J162 1, In each panel the black spectrum represents the lensed 
galaxy and the composite spectrum of 81 LBGs from the analysis of Jones ct al. (2012) is shown in red. All spectra are plotted in the rest frame 
with flux normalized such that median /„ = 1 in the range 1250 — 1500 A. Prominent spectral features are labeled. 
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Fig. 2. — Mean low-ionization absorption line profile and its associated neutral gas covering fraction derived using the methods discussed in 
Section [3] Si II covering fraction of the z ~ 4 arcs is measured from smoothed spectra with FWHM resolution ~ 110 kms , while average 
line profiles are from unsmoothed data (FWHM ~ 70 kms -1 ). There is no si gnific ant difference in results derived from the smoothed and 
unsmoothed spectra. Velocities are relative to adopted systemic redshifts (Scction l3,l| l. derived from absorption line centroids for A2390-H3 and 
nebular emission in the othe r two cases. Equiva lent measurements from ESI sp ectra of z = 2—3 arcs are shown for comparison. Details of these 
ESI spectra can be found in lPettini e~ al. (2002) and Quid er et aTl ((2009, 2010). 
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Fig. 3. — Outflow kinematics and covering fraction as a function of redshift. Top: maximum blueshifted velocity at which outflows are 
detected. In the case of A2390-H3 there is an uncertainty of ~ 125 kms -1 arising from the unknown systemic redshift. Bottom: maximum depth 
of the average absorption profiles shown in Figure [2] whic h serves as a proxy for the covering fraction of H I at the corresponding velocity. We 
include estimated v alues for the 8 o'clock arc (Dcssaugcs-Zavadsky ct al. 201(]) and a somewhat lower resolution (R=1300) spectrum of Q0000-D6 
USteidel ct al. 2010) in addition to the sources shown in FigurefJ] Galaxies are color-coded according to their H^Lya and relevant data are listed 
in Table [2] On average, the z = 4 galaxies show somewhat lower maximum outflow velocities and lower covering fractions (weaker absorption 
depth) compared to similarly studied sources at lower redshift, but the strongest trend is a decreasing covering fraction with WL ya . The 2 = 4 
data enable us to exa mine trends of H I covering fraction at a time significantly closer to the epoch of reionization, thought to end at z ~ 7 (e.g., 
ISchenk cr ct al. 20lJ). 
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Fig. 4. — Correlation of Lya equivalent width with H I covering fraction as traced by low-ionization absorption lines. The dashed line shows a 
linear fit to the data. Here WLya includes only the emission component, and the two points at (1,1) are slightly offset for clarity. As discussed in 
the text, thes e absorption line measur ements correspond to a lower limit on f c and an upper limit on f eS c- Since Lya emission strength increases 
with redshift (Stark ct al. 2010, 2011), this result likely indicates lower covering fractions (which would permit higher escape fractions) at earlier 
times. 



